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Abstract 

We prove that in random matrix theory there exists a universal 
relation between the one-point Green's function G and the connected 
two- point Green's function G c given by 

N 2 G c (z,w) = Qzdw ^°s(( G ^zZ^"'^ ) + irrelevant factorized terms. 
This relation is universal in the sense that it does not depend on the 
probability distribution of the random matrices for a broad class of 
distributions, even though G is known to depend on the probability 
distribution in detail. The universality discussed here represents a 
different statement than the universality we discovered a couple of 



years ago, which states that a 2 G c (az, aw) is independent of the prob- 
ability distribution, where a denotes the width of the spectrum and 
depends sensitively on the probability distribution. It is shown that 
the universality proved here also holds for the more general problem 
of a Hamiltonian consisting of the sum of a deterministic term and a 
random term analyzed perturbatively by Brezin, Hikami, and Zee. 
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1 Introduction 



The theory of large random matrices |5|, |6], [7|, [27[] has been extensively de- 
veloped over the years. Recently, in a series of papers [|], || 0, [| |30l1 , 
we, and with the collaboration of J. D'Anna and of S. Hikami, have studied 
the correlation between the density of eigenvalues of large random matrices. 
For physical applications, one may imagine under some circumstances repre- 
senting the Hamiltonian of a disordered system by a large hermitean random 
matrix. 

We will try to make this paper as self-contained as possible and so we 
will begin with the necessary definitions. Denote by ip an N by N hermitean 
matrix taken from the probability distribution 

p(^) = }_ e -mrv{<p) (1) 

Zj 

with Z fixed by / d<pP(tp) = 1. Here V(<p) is an arbitrary polynomial of tp 
and the large N limit is understood. A hermitean matrix <p is picked with 
the probability P((p) and we imagine that its eigenvalues are found. When 
this procedure is repeated a large number of times, we can define a density 
of eigenvalues, averaged over the matrices in an ensemble defined by the 
distribution ([[]). Similarly, the density- density correlation function can also 
be defined and calculated. 

Before we continue with the necessary definitions, let us explain the phi- 
losophy underlying our work. We have devoted considerable effort to dis- 
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covering and elucidating universality properties satisfied by the correlation 
function, and so it behooves us to say a few words about our motivations. 
We say a quantity is universal if it does not depend on the probability dis- 
tribution P, or for the specific class of P examined here, on V. 

While there are notable exceptions of course, much of the literature on 
random matrix theory we are aware of is concerned with only the Gaussian 
distribution, namely when V(ip) is quadratic in ip. While this assumption 
may be plausible, it is important to understand to what extent the results 
obtained are general. Given a Gaussian distribution, essentially everything 
can be calculated. In particular, in the orthogonal polynomial approach, the 
relevant polynomials are just the Hermite polynomials. When asked why 
they restrict themselves to the Gaussian distribution, many workers in this 
subject would simply assert that it "suffices." Others take the attitude that 
one does what is possible, and indeed if one were to go beyond random 
matrix theory to study actual disordered systems, it is difficult to treat non- 
Gaussian disorder. 

What we have been trying to do is to see, at least within the confines of 
random matrix theory, which results are indeed universal and which results 
are specific to the Gaussian case. This search for universal quantities is made 
all the more important by the fact that the density of eigenvalues is in fact 
known to be not universal in the sense used here, that is, the density of 
eigenvalues depends sensitively on V, as will be mentioned below. 
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We would like to distinguish the universality discussed here from the 
"short distance" universality quoted in the literature. By short distance, we 
mean that we look at the correlation function on scales larger than, but com- 
parable to, the spacing between the eigenvalues. On this scale, the physics is 
essentially controlled by level repulsion, leading to a more or less locally uni- 
form spacing between eigenvalues. Thus, we expect the correlation function 
to be universal when suitably scaled by the local density (see for example 
equation (2.19) of [0). The universality which we have discussed in our work 
and which we will discuss here is on distance scale large compared with the 
spacing between the eigenvalues and is thus less evident. 

Universality in critical phenomenon is of course well understood by now, 
essentially because of the diverging correlation length near a second order 
phase transition. The subject has been placed on a solid foundation by the 
renormalization group approach. In [|I|] we suggested that a renormalization 
group inspired argument may also be operative here, although the details 
are still not completely clear, particularly since the density of eigenvalues is 
not universal. Perhaps, using the language of renormalization group, we may 
suppose that a "relevant operator" is involved. 

Now that we have stated some of our motivations let us continue with 
our definitions. We introduce the Green's function (or resolvent) 

The bracket will henceforth denote averaging with respect to P: for any 0(ip) 
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we write < 0((p) >= J dipP(<p)0(<p). 

The density of eigenvalues of the random matrix <p is then given by p(p) = 
-^tr<5(p — <p)\ = — ^ImG(p + ie). The limit N tending to infinity is always 
understood. Notice that in this paper, as in our earlier work, we choose the 
factors of N in our various definitions such that the interval over which p(p) 
is non-zero is finite (i.e., of order iV ) in the large N limit. 
The two-point Green's function is defined by 

G(z, w) = (hr—Lr—) (3) 
\JM z — (f JM w — ip I 

In the large limit, G(z, w) — > G(z)G(w) and thus it is customary to define 
the connected Green's function defined by G c (z,w) = G(z,w) — G(z)G(w), 
a quantity of order 1/N 2 . (Thus, we will be dealing with a quantity often 
ignored in discussions of large N expansion.) The connected correlation 
between the density of eigenvalues is then given by 



c 



= (-1/4tt 2 )(G c (++) + G c (— ) - G c (+-) - G c (-+)) (4) 

with the obvious notation G c (±, ±) = G c (p±ie, v±i5) (signs uncorrelated) . 
Almost twenty years ago Brezin, Itzykson, Parisi, and Zuber [f^| calcu- 



lated the one-point Green's function and found that, as might be expected, it 
depends on V in a complicated way. Purely for the sake of completeness, let 
us record that for V(<p) = J2l=i lk~9 ki P 2k ( we ^ a ^ e ^ ^° ^ e an even polynomial 
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for simplicity) the Green's function has the form 



G(z) = -\V\z) - P(z)Vz^^} (5) 



where the polynomial 



p^ = ltmt { M4r^- 2n - 2 (6) 

Z k=l n=0 V b -) * 



The endpoint a of the spectrum is determined by 



If a (2A0! f «V-l 



(7) 



The density of eigenvalues is given by 



p(/x) = -P(/x)^a 2 -/i 2 (8) 

In particular, in the simplest case, with the Gaussian distribution defined by 
V(ip) = m 2 (p 2 /2, the one-point Green's function is given by 



Go(z) = -[z - V#=Z] (9) 

and thus the density obeys Wigner's celebrated semi-circle law 

M = A^a 2 -^. (10) 

We will emphatically not need any of these explicit formulas in what 
follows. We simply want to emphasize to the reader that, not surprisingly, 
G(z) and p(/x) both depend on the potential V in detail. We say that the 
one-point Green's function and the density are not universal. 



In [|T| we calculated the correlation function p c (/i, v) using the method 
of orthogonal polynomials and again as might be expected found that it 
depended in detail on V. We will not display the full expression here but 
simply note that the resulting p c (/i, v) oscillates wildly as a function of \i and 
v on scales of 1/N. These oscillations are entirely expected since between 
/i and v finitely separated (that is, with // — v ~ 0(1)) there are in general 
O(N) eigenvalues. Thus, it is natural to smooth p c (/i, v) by integrating over 
intervals b\i and 5u, large compared to 0(A^ _1 ) but small compared to 0(N°), 
centered around fi and v respectively. Upon smoothing, the full expression 
for the correlation function simplified enormously and we obtained |l| 

smooth/ x = _Z1 I ~ (11) 

Pc ^' J 2N 2 n 2 {n-v) 2 [{a 2 - fi 2 ){a 2 -v 2 )} 1 / 2 ' 1 ' 

We found that, remarkably enough, the smoothed correlation function de- 
pended on the polynomial V only through the single quantity a, the width of 
the spectrum. In other words, if we introduce the obvious scaling variables 
x = fx/ a and y = vja then the correlation function (henceforth we will drop 
the superscript "smooth") is equal to 



a 



with the universal function 



f( \ 1 (1 ~ xy) 

J{X,y> {x - y) 2 [(1 - x 2 ){\ - y 2 )f/ 2 1 ' 

This universality has since been derived by other authors |2^, using 
alternative methods, and verified numerically plj. It is perhaps useful to 



remark here that, while there are notable exceptions of course, much of the 
literature on random matrix theory, as far as we know, is devoted to the 
Gaussian case |12, 13| . The whole point of our work is that it is possible to 



go beyond the Gaussian distribution. 

In [[§ we developed a diagrammatic approach to calculating the connected 
two-point Green's function G c (z, w). We will describe this diagrammatic 
approach in detail in the next section. Here we will simply outline our result 
from ||. (To read the rest of this paper, it is not necessary to have read 
[||] first.) Using a diagrammatic approach, we find that G c (z, w) is given by 
an infinite set of Feynman diagrams. We were able to calculate and sum 
this infinite set only for the Gaussian case and obtained a relatively simple 
expression for G c (z,w) (see equation (2.11) in [[|). Later we recognized[|] 
that this expression may be written in the elegantly compact form 

iV 2 G.^, m ) = ^l°g( GoW " GoW ) (14) 
ozow z — w 

The subscript "0" indicates that the quantities in this equation are all cal- 
culated with the Gaussian distribution. Taking the absorptive part of this 
according to (Q) we obtain the smoothed correlation function given in ([12]). 
We explained in 0] that the diagrammatic method "automatically" gives the 
smoothed correlation function. This is because in the diagrammatic method 
we calculate the Green's function G c (z, w), by first letting N go to infinity to 
pick out an appropriate set of diagrams, and then letting z and w approach 
the real axis to extract the correlation function from the absorptive part of 
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G c (z, w) according to (f|). With N going to infinity, the discrete set of poles 
of G c (z, w) on the real axis merges into a cut. This is equivalent to the 
smoothing procedure employed in 0]. 

In ||, in contrast to the work we did in 0, we were unable to calculate 
G c (z, w) for a general V. Indeed, the task of summing up the infinite sets of 
graphs generated by the interaction terms in V appeared to us at the time 
enormously complicated and perhaps even hopeless. It is simple enough, 
however, to summarize the remarkable universality discovered in [|TJ. The 
universality expressed in fll2|) and fll3|) can be stated in terms of G c (z,w) by 
saying that for a general V we have 

^C efero ) = ^-log( G ° (z/a) - G ° (m/a) ) (15) 
ozow z — w 

It is far from obvious how such a relation can be derived diagrammatically. 
Here a has a complicated dependence on V as indicated by (|7|). Furthermore, 
the Gaussian Green's function Gq(z) appears on the right hand side. The 
Gaussian Green's function Gq(z), in contrast to the Green's function G{z) 
appropriate to the general V, is not a "natural" object to appear in a calcu- 
lation of G c (z,w). Yet, according to the orthogonal polynomial analysis of 
[U], this relation must be true! 

To appreciate how complicated a diagrammatic calculation of G c (z,w) 
can get, the reader is invited to look at || where together with Hikami we 
attempted this calculation. We had to restrict ourselves to the gip 4 case, and 
even so, we were able to obtain, after a long and rather involved calculation, 
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the correlation function G c (z,w) to only first order in g. This calculation, 
however, was instructive. It turned out that the numerous terms in our fi- 
nal expression for the two-point Green's function G c (z, w) can be grouped 
together in precisely such a way that the rather complicated final expression 
can be written in terms of the one-point Green's function G(z). This ob- 
servation is highly non-trivial in that, as we can see from (||), @, (]?]), the 
Green's function G(z) to first order in g is already not particularly simple. 
Brezin, Hikami, and Zee|| found that, to first order in g, 

N*G c (z, w) = JfL \og(( G{z) ~ G{W) ) (1 + AgG(z)G(w))- 1 ) + 0{g 2 ) (16) 
ozow z — w 

Note that the factor of (1 + 4gG(z)G(w)) to this order in g contributes to 
G c (z, w) only a factorized term like h(z)h(w) for some function h. These 
factorized terms would not contribute to the connected correlation function 

It was thus tempting to conjecture that for an arbitrary V(<p) the con- 
nected two-point Green's function can be written as 

n , , 9^ . ,G(z \ — Giw\ , . , , 

JM G c (z,w) = logff ) + irrelevant factorized terms (17) 

ozow z — w 



For the Gaussian case, the two expressions in (|15|) and ([17|) are manifestly 
the same. For a general V, however, the equality of these two expressions is 
far from evident. 

In a recent paper 0] we adopted a slightly different philosophy: instead 
of calculating G(z) and G c (z, w) in terms of V and then trying to express 
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G c (z, w) in terms of G(z) and G(w) by eliminating the dependence on V, we 
attempted to calculate G c (z, w) directly in terms of G(z) and G(w), appealing 
to V only for the general structure of the Feynman diagrams. This shift in 
philosophy is reminiscent of the dispersion approach in particle physics in the 
1950's: instead of trying to calculate various physical quantities in terms of 
an underlying Lagrangian, particle physicists of that era attempted to relate 
various physical quantities to each other. In this paper, we will exploit this 
philosophy to prove the conjecture in (jl7f) . In line with this philosophy, we 
will keep the amount of explicit calculation to a minimum. Instead, we will 
organize the relevant diagrams in such a way as to obtain structural relations 
between different Green's functions. 

2 Diagrams 

Let us review the diagrammatic approach discussed in || . We may regard the 
distribution ([!]) as defining a (0 + 0)-dimensional field theory. In this context 
the Feynman diagram approach consists of nothing more than expanding 
G(z) in inverse powers of z and doing the integrals in 

OO 1 1 

G(z) = E < jjW > (18) 

n=0 A 1 V 

In doing the integral over (p implied by < > we split off the quadratic 

part V(<p) and treat the rest of V perturbatively. As explained in 0, it is 
useful to borrow the terminology of large N quantum chromodynamics 



from the particle physics literature, and speak of quark and gluon lines. See 
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figure (1) for a graphical representation. (It is of course not necessary to use 
this language, and readers not familiar with this language can simply think of 
the diagrams as representing the different terms one encounters in doing the 
integral in @.) The bare quark propagator simply comes from the explicit 
factors of z in ([0|) and is represented by a single line and given by 1/z. The 
quadratic term in V(ip) determines the bare gluon propagator, represented 
by double lines, and proportional to 

(*>>?) = ti^ (19) 

where m 2 is defined by the quadratic part of V((p) = m 2 ip 2 /2 + The 

non-Gaussian terms in V(ip) describe the interactions between gluons. 



The important point, as originally stressed by 't Hooft [23], is that this 
double-line formalism provides an efficient way of counting the powers of N 
to be associated with each diagram: each vertex counts for one power of N, 
each gluon propagators counts for iV -1 , and each closed loop counts for N. 

The bare quark propagator 1/z is changed by the interaction to the 
dressed quark propagator G(z). The gluon propagator is dressed by gluon 
interaction, but note that it is not dressed by quark loops. This is clear from 
the definition of our problem. Another way of saying this is to note that the 
one-point Green's function may be represented, by using the replica trick, as 

G(z) = Um n ^o I D^Di/jDipP^ipl^e-^i^-^ (20) 



The replica index a runs from 1 to n. Note that in this language the ^'s 
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represent the quark fields and ip the gluon fields. The interaction between 
gluon and quarks are given by ipaftpa- (Color indices are suppressed here.) 
The interaction of the gluons with each other is determined by the non- 
Gaussian part of P(<p)- Since internal quark loops are proportional to the 
number of replicas n, they vanish in the n — > limit. 

Similarly, we can treat the two point Green's function by expanding 

G ^ H S£i^^(^V^") (21) 

The implied integration over ip then generates the Feynman diagrams for 
G c (z, w). We can also use the replica trick to represent G c (z, w). Clearly, we 
would have to introduce two quark fields ip z and ip w : the variables z and w 
act like a flavor label. Thus, the correlation function describes two quarks, 
"carrying" z and w respectively, interacting by emitting and absorbing glu- 
ons (which have complicated interactions amongst themselves.) What we 
are doing here may be considered as a "baby version" of quantum chromo- 
dynamics. 

As a warm up exercise and to gain some familiarity with what is going 
on, we will first consider the Gaussian case. The derivation given here is 
simpler than the one given in || and in its essence was given in one of our 



earlier papers | 15| . With the benefit of hindsight, we start by taking out two 



partial derivatives: 



. / 1 11 1 
& c (z,w) = ( — tr —tr- 



N z — ip N w — ip ! 
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= d- Z d^(w Mog{z - ^ trl0 ^ W - ( 22 ) 
Expanding the logarithms, we find 

This is represented by the "wheel" graph of Fig (2) where the quark propa- 
gator on the inner rim carries z and the one on the outer rim carries w. For 
the moment we ignore quark self-energy and vertex corrections: every gluon 
emitted on the inner rim is absorbed on the outer rim, and vice versa. Thus 
we may set n = k in (0). 

Since we are working with a Gaussian distribution we can immediately 
evaluate (tr(p n tr(p n ) c = n. (Here with no loss of generality we have scaled 
m 2 to unity.) Graphically the factor of n corresponds to the fact that with 
the inner rim held fixed, we may rotate the outer rim by n different "clicks" 
and leave the diagram invariant. It is this factor of n which produces the 
logarithmic function when we evaluate the sum in ( p3| ) to obtain 

N 2 G(z, w) c = ~ A log(l - — ) (24) 
oz aw zw 

This expression does not yet have the form in flT7D. Next we have to include 
self-energy and vertex corrections. Instead of doing this let us leave this 
expression as it is for the moment and turn our attention to a "scattering" 
formalism discussed in [ [TEH . 
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3 Some formalism 



Let us go back to the expansion of G c (z,w) in (|2lD . Note that there is an 
extra power of 1/z and 1/w compared to the powers of (p. Thus, in the 
wheel diagram of figure (2), one of the quark propagators on the inner rim, 
and one on the outer rim, should actually be represented by 1/z 2 and 1/w 2 
respectively: they each consists of two quark propagators. We represent this 
fact graphically by two dots, one on the inner rim, and one on the outer rim, 
as shown in figure (3). There are two possibilities: the two dots are on the 
same "sector" of the wheel, as shown in figure (3a), or the two dots are on 
different "sectors", as shown in figure (3b). 

Now imagine cutting open the quark propagators at the two dots. This 
converts the wheel diagrams into two sets of scattering diagrams as shown in 
figure (4). Note that it is necessary to include the crossed ladders in figure 
(4b). (Incidentally, the necessity of including the crossed ladders came to us 
as a bit of a surprise in carrying out this calculation using the formalism of 
U but it is made completely clear by the present formalism.) 

4 Some formal relationships 

We have illustrated the discussion in the two preceding paragraphs with 
diagrams appropriate to the Gaussian case, but this discussion applies im- 
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mediately to the general case. Let us define the "scattering amplitude" 

n((— )){-!— )Z ) = N(((—))(—)A - 5)G{z)8™G{w)) 
\ z — <p J w — <p I \ z — ip J w — (p I J 

= + 8}8™B (25) 



The two "scalar" scattering amplitudes A and B depend on z and w of 
course and correspond diagrammatically to the sets of graphs shown in figure 
(5a, b, c). We see from the flow of the color indices that A and B correspond 
to the two ways of cutting the wheel diagram, that is, to figure (4a) and 
figure (4b) respectively 

For a general V it is complicated to calculate the amplitudes A and B 
directly by perturbation theory. For example, A is given by the infinite set 
of graphs in figure (5a) for the gf A theory. However, we will see that by 
a judicious arrangement of our calculation, we can avoid doing the explicit 
calculation that we had to work hard to carry out in simply to obtain a 
result to first order in g. 

First, we note that by contracting (25) with 5^5^ we have 

N 2 A + B = N 2 [ G{Z) ~_ G{W) - G(z)G(w)] (26) 

We can check easily by looking at a few graphs that A is of order iV° while 
B is of order 1/N. Thus, we can drop B in this equation and determine A 
(to leading order in N of course) in terms of the one-point Green's function 
G(z). 
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Next, by contracting (|25| ) with Sjd 1 ^ we find that the connected two-point 
Green's function is given by 

N 2 G c (z, w) = ltr—^—tr——) =A + NB (27) 
\ z-(f w-y/ c 

Note that in this equation A and B both contributes to the same order in 
N. Thus, to determine G c we still have to know A and B. It would seem 
that we would have to work to obtain B, as we did in our previous papers. 
In the next section, we will see how we can avoid calculating B. 

5 General V 

We are now ready to tackle the full problem of determining G c (z, w) for a 
general V. First, it is useful to define a two-quark irreducible scattering 
amplitude T(z, w) consisting of those graphs that do not fall into two discon- 
nected pieces upon cutting the two separate quark propagators, as shown in 
figure (6). Then the scattering function A is evidently given by 

a = J— (r + r— r + r— r— r + ) = _s£_ (28) 

z z w z zw zw zw 1 1 

zw 

Note that this expression merely relates A to T(z,w), which at this stage 
would appear to be an exceedingly complicated object to calculate directly 
Let us now start the computation of G c (z,w) with (^l|) which we repeat 
here for convenience: 

OO OO 1 , -1 1 , 



Again, let us proceed by first ignoring vertex and self energy corrections. But 
in the general case we can no longer simply set n equal torn, as we did in the 
Gaussian case. We see, however, that we can express all the wheel diagrams 
representing (|2ll) in terms of the (unknown) amplitude T(z,w) as indicated 
in figure (7). 

We have yet to put on the two dots, one on the inner rim, one on the outer 
rim, as explained above. We see that now there are a number of possibilities. 
We can put the dot on an "exposed" quark line, or we can put the dot on a 
quark line hidden inside a T. These two possibilities are illustrated in figure 
(8). 

For obvious reasons, we now find it useful to use an alternative notation 
in which we replace z and w by z x and z 2 respectively and define d a =-^-. 
When we put the dot on an exposed quark line, we replace l/z a by l/z a 2 = 
—d a (l/z a ). When we put the dot on a "hidden" quark line, we in effect 
replaced T(zi,z 2 ) by — d a T(zi, z 2 ). Incidentally, in the Gaussian case there 
is no "hidden" quark line: all quark lines are exposed by definition. The 
differential operator <9i<9 2 associated with putting on the dots is precisely 
what relates flU) to (||). 

In addition to the choice of putting the two dots on exposed or hidden 
quark lines, we also have the choice of putting the two dots in the same 
sector or in different sectors. (Sectors are defined as the segments of the 
wheel divided by the different T's: each sector consists of one the spaces 
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between the T's and a T next to that space.) 

Let us first consider putting the two dots in the same sector. Then we 
have the following four possibilities, corresponding to the four diagrams in 
figure (9): (a) both dots are on an exposed line, thus giving -f-y, (b) the dot 
on the inner rim is on a hidden line, while the dot on the outer rim is on 
an exposed line, thus giving = ^-, (c) the previous case with inner and outer 
exchanged, and (d) both dots are on a hidden line, thus giving 8 * ^ r . These 
four terms add up to + =g£ + ^ + = ^(^). 

The rest of the wheel (see figure (9)) can be filled with nothing, one T, 

two T's, and so on, that is, with the series 1 H — - — h (-^— ) 2 + Putting 

all of this together we have the following contribution to G c (z 1: z 2 ): 

{—^nd^—) (so) 

J- ~ TT~ z l z 2 

Z\ Z2 

Next we have to consider the possibilities of putting the two dots on two 
different sectors. The dot on the inner rim can be either on an exposed line or 
a hidden line, and thus we obtain a factor (-£ — = —di(-^—). Similarly, 

' V 2^22 2122 ' ±v 2l22 / J 1 

the dot on the outer rim can be either on an exposed line or a hidden line, 
and we obtain the factor just given but with 1 and 2 interchanged. These 
two factors combine to give di(-^—)d 2 (-^) 

The two different sectors, where the two dots are placed, divide the wheel 
into two segments, each of which can be filled, just as above, with nothing, 
one T, two r's, and so on, that is, each segment leads to the factor - \ . 

z l z 2 
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Thus, we obtain in G c (z\,Z2) the contribution 



2 d x {^-)d 2 {^-) (31) 



i' ' ZiZ2 ' ZiZ2 



Putting ( j30|) and fl3l|) together, we finally obtain 

1 , r s 



N 2 G c (z u z 2 ) = (- H^ 2 ( 



1 "1717 ^1^2 ' 



-122 



which we happily recognize as just 

N 2 G c (z u 22) = -^1^2 log(l - — ) (33) 

-21-22 

Note that we could have given a shorter derivation by simply "working 
backwards": we could have started with (^) and simply said that the op- 
erator (dialog) distributes the two dots in all the possible ways that we 
had enumerated, but we believe that our longer derivation just given is more 
transparent and easier for the reader to follow. 



For T = 1 we recover our previous expression (24j) : surely we are on 
the right track. We are almost there but we have yet to put in the quark 
self energy corrections and the vertex corrections. The quark self energy 
corrections are easy to put in: we simply replace the bare quark propagator 
1/z by the dressed propagator G(z) appropriate for the interaction potential 
V. The vertex corrections require more thought. First, note that our usage 
of the term "vertex corrections" differs slightly from the standard usage. For 
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example, the diagram in figure (10a) has already been counted in T. We 
include in vertex corrections the diagrams in figure (10b) for example. As 
seen in figure (10c), including the vertex corrections we simply multiply the 
amplitude without vertex corrections by a factor v(z). 

Already, we noted in || that the gluons in the vertex corrections must 
"span" the whole amplitude |20 lest we lose factors of N. This is illustrated 



in figure (11). The remark given in || concerning the vertex corrections for 
the Gaussian case clearly generalizes to the case of an arbitrary V. 

We are now faced with the task of calculating v(z), which we calculated 
for the Gaussian case in ||. Fortunately, we can avoid doing any work by 
noting that there is a 'Ward identity" 

„ W = ^ (34) 

The reader can easily convince himself or herself of this identity by contem- 
plating the diagrams in figure (12). The differentiation ^ simply puts the 
vertex in all possible places with the correct counting factor. Note that the 
v(z) given by this identity is fully dressed, that is, the quark propagators 
that enter in v(z) are already dressed. 

But now we see a remarkable cancellation of the vertex corrections be- 
cause we have the foresight (or hindsight!) of arranging our calculation of 
G c so that it has the form d\ 82 {■■■■) before self energy and vertex corrections 
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are included (see (p^)). Including these corrections we obtain 

r „, 



N'G c (z x ,z 2 ) = v(z 1 )v(z 2 )([d l d 2 \og(l 



Z\Z 2 



\dressed) 



r 

= <9i<9 2 (log(l \dressed}) (35) 

Z\Z 2 

since {£)\dressed = dG -i( z ) = ^ys- Tne vertex corrections disappear! 



At this point, if we go back to the Gaussian formula ( ^ip we see that in 
the Gaussian case we have essentially finished our calculation since T = 1. 
We obtain 

N 2 G Qc (z, w) = -—— log(l - G Q (z)G (w)) (36) 

Inserting the explicit form for Gq(z) given in (|9]) we obtain after some simple 
manipulations (P^). This derivation is simpler than that given in (0). 

It is always nice to recover the Gaussian result as a check but here we 
want to do the much more ambitious problem of calculating G c for a general 
V. To go further, we have to calculate T(zi,z 2 ) and then to dress it by 
replacing l/z a by G(z a ) which we will write as G a for short. This would 
have been a long involved calculation, but again we note happily that we can 
avoid doing it simply by noting that the bare T is related to the bare A by 



r = -^4- (37) 



But the dressed A is given by the identity in (|^)! Thus, we don't have to 
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do any further work. We simply dress fl3~T|) to obtain 



^dressed ~ ( ^ n H ) (38) 



Finally, then 



[Z\Z2 — Z2))dressed — ( n n ^dressed) — ~ (7; T=r) (39) 



And thus, we have proved our conjecture. Combining (|35| ) and ( |39D we obtain 
our conjectured relation 

N 2 GJz, w) = - - log(( ) + irrelevant factorized terms (40) 

ozow z — w 

Finally, in the language of the wheel diagram, it is easy to see where the 

irrelevant factorized terms in ( HOI) come from. They come from diagrams 

which disconnect the inner rim and the outer rim of the wheel from each 

other and thus clearly has a factorized dependence on z and w. See figure 

(13). 

6 Deterministic plus random 

In our earlier work, we have also generalized the problem outlined in the 
introduction of this paper to the problem of a Hamiltonian given by the sum 
of a deterministic term and a random term 

H = H + v (41) 

Here H is a diagonal matrix with diagonal elements % = 1, 2, ...N, and (p a 
random matrix taken from the ensemble ([]]). For the Gaussian case, namely 
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with V(<p) = ^(p 2 , Pastur [^J has long ago determined the density of eigen- 
values. Our work in [R) went beyond Pastur's work in that the correlation 
function between the density of eigenvalues in the Gaussian case was also 
determined. In our recent work with Hikami, |§ we studied this correlation 
function for a g(p 4 theory to first order in g. 

This problem of "determinism plus chance" may be regarded as a generic 
problem in physics, and as such represents a significant generalization of 
Wigner's problem. For example, consider an electron moving in a magnetic 
field and scattering off impurities. We note that these "deterministic plus 
random" problems are considerably more difficult than purely random prob- 
lems. The orthogonal polynomial approach used in [|IJ involves diagonalizing 
the random matrix ip and is clearly no longer available: in ([0]) we cannot 



diagonalize ip without un-diagonalizing Hq. Thus, we do not have the analog 
of flT5|) for this problem. 



Indeed, as discussed in a recent paper J8|, the problem described here rep- 
resents a special case of a broader class of problems involving the addition of 
random matrices. The deterministic Hamiltonian H Q may in turn be replaced 
by a random Hamiltonian. Indeed, a deterministic matrix is but a special 
case of a random matrix. 

Consider then a Hamiltonian given by 

H = <p! + <p 2 (42) 
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with the matrices tpi^ taken from a factorized probability distribution 

P(<Pi, V2) = le-^IViC^O+v,^)] _ p l(v9l )p 2 (^). (43) 

In [U it was shown how the one-point Green's function G can be obtained 
for the Hamiltonian given in fl42|) . Here we would like to solve the problem 
of determing the connected two-point Green's function G c . 

A slightly sloppy but essentially correct argument is that given our uni- 
versal relation (^) between G c and G our problem is solved instantly. The 
desired connected two-point Green's function is given in terms of the one- 
point Green's function appropriate to the distribution in (|43|). This result is 
precisely what was conjectured in ||. 

We can put this argument on a more solid footing by using the formalism 
discussed in M. (The following discussion will be sketchy and not self con- 
tained.) In that work, it was shown that G(z) may be determined in terms 
of G\(z) and G 2 (z), the Green's functions corresponding to the distribution 
Pi and P 2 respectively, according to the following procedure. First, solve the 
equations G a (B a (z)) = z, for a = 1, 2, that is, find the functional inverses of 
G a (z), denoted by B a (z) here. Next, define the function 

B 1+2 (z) = B 1 (z)+B 2 (z)-\ (44) 

The functional inverse of -81+2(2;) is then the desired Green's function G(z). 
This type of addition laws has been discussed recently in the mathematical 
H], [TJ|, [L9| and physical literature [|14], [LJ], [L6| . 
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Following the argument given in M which we won't repeat here, we find 
for example the undressed T is given by 

n/ \ /-, z 2 G gcl (z) - w 2 G gcl (w) 

T(z, w) = zw(l ) + (1 <-» 2) (45) 

Here G gc \ and G gC 2 are "gluon connected" Green's functions defined in || 
and are related to B\ and B 2 respectively. Following the same steps as above 
we find that fl35|) still holds. Inserting the expression for V given here and 
using ([H]) we immediately find that ( |40D indeed holds for this more general 
class of problems. 

7 Conclusion 

In conclusion, we have found a remarkable universal relation between the one 
point Green's function G(z) and the connected two point Green's function 
G c (z,w). This represents an entirely different sort of universality as the one 
found in There it was shown that the scaled two point Green's function 
a 2 G c (az, aw), with a the endpoint of the spectrum given by a complicated 
function of the potential V(ip), is independent of V. Here it is shown that 
the structural relation between G c (z,w) and G(z) is independent of V, even 
though G(z) is known to depend on V in a complicated way. 

While we know that these two forms of universality must be equivalent, 
it is not obvious how to show this equivalence directly. 

The compact form of G c (z, w) obtained here renders the universality prop- 
erty of p s c mooth (iJ, v) as approaches v particularly transparent. Consider our 
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universal form 

Ir9 ^ , . d d . ,G(z) — G(w). ,,„ s 

N 2 G c {z, w) = — — log ^ 46 

ozow z — w 

If z and w approach each other on the same side of the cut of G, the argument 
of the logarithm is a smooth function of z — w, and thus would not contribute 
to p s c mooth (n, v) a term proportional to — v) 2 that we know from (|Tl"l) 
must be there. On the other hand, if z and w approach each other from 
opposite sides of the cut, then writing G(fx ± ie) = R(fJ,) ± il(fj) as e goes to 
zero, we have the universal singular behavior 

N 2 G c {fi + ie,v-i8) = — — log(^^ K -L + i^- 

9 d . . . 1 

Inserting this into (ffl) we find the singular part of ( ]TT| ) . The emergence of the 
universal behavior, in which the dependence on R(fi) and I(^) drops away, 
is made particularly clear by ([47]). 

We may be tempted to conjecture that the elegantly compact form of 
our universal relation (f46|) may be associated with a deeper mathematical 
structure. 
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9 Figure Captions 

Fig 1. Diagrammatic rules: (a) quark propagator, (b) gluon propagator, (c) 
quark gluon vertex, and (d) gluon interaction, illustrated here with a gtp A 
vertex. 

Fig 2. Wheel diagram. 

Fig 3. The two dots may be placed in the same sector (a) or in different 
sectors (b). 

Fig 4. The diagrams in Fig 3. cut open at the dots. (The diagrams in the 
same topological class, rather than the exact correspondents, are shown.) 

Fig 5. Diagrams contributing to A (a) and diagrams contributing to B 
(b) and (c). Note that two topological distinct classes of graphs contribute 
to B. 

Fig 6. Typical diagrams contributing to T. 

Fig 7. The wheel diagrams for a general potential V. 

Fig 8. A dot may be put on an exposed line (a) or on a hidden line (b). 

Fig 9. (a) The two dots are both placed on exposed lines, (b) The dot on 
the outer rim is placed on an exposed line, while the dot on the inner rim is 
placed on a hidden line, (c) The situation in (b) reversed, (d) The two dots 
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are both placed on hidden lines. 

Fig 10. (a) Diagram not included in what we called vertex correction, (b) 
Diagram included in what we called vertex correction, (c) A vertex correc- 
tion: the shaded portion can include many gluon lines, possibly interacting 
with each other. 

Fig 11. (a) The gluon in the vertex correction spans the whole diagram, 
(b) The gluon in the vertex correction spans only part of the diagram. We 
see that the diagram in (b) has one less loop than the one in (a). 

Fig 12. The "Ward identity" relating the vertex correction to the self- 
energy. 

Fig 13. A diagram contributing to the "irrelevant" factorized terms in 
N 2 G c (z, w). Note the inner rim and the outer rim are "decoupled" from each 
other. 
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